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Abstract 



A regularization for the baryon number consistent with the energy in the Nambu-Jona- 
Lasinio model is introduced. The soliton solution is constructed with the regularized baryon 
number constrained to unity. It is furthermore demonstrated that this constraint prevents 
the soliton from collapsing when scalar fields are allowed to be space dependent. In this 
scheme the scalar fields actually vanish at the origin reflecting a partial restoration of chiral 
symmetry. Also the influence of this constraint on some static properties of baryons is 
discussed. 
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Introduction 



In recent years the investigation of the solitons of the Nambu-Jona-Lasinio (NJL) model 
[|l| have experienced steady progress 0, [3|, |J ||. While in the Skyrme model |7|, |8j the 
baryon number is assumed to be given by the topological charge the chiral soliton of the 
NJL model is not strictly a topological soliton. On the contrary, the baryon number is given 
more directly in terms of the quark fields and, in general, contains contributions from the 
valence quarks as well as from the Dirac sea . In leading order of the derivative expansion 
the latter is given by the topological charge ||11|| . Whether the baryon number is dominantly 
carried by the sea or the valence quarks depends on the specific features of the model ||. 

Since the NJL action is not renormalizable regularization is required and more impor- 
tantly the resulting energy functional, baryon number density, charge radii, etc. will depend 
on the regularization scheme employed |12[. Regularization is commonly applied to the 
Wick-rotated Euclidean NJL action. A special feature of this Euclidean action is the fact 
that only its real part is ultraviolet divergent while the imaginary part stays finite. Ac- 
cordingly in most treatments only the real part undergoes regularization. As long as time 
components of vector and axialvector fields are ignored only the real part contributes to the 
static energy functional, falsely pretending that the imaginary part is of no importance for 
soliton solutions. The imaginary part is indeed relevant because it completely determines 
the baryon number current. Since we intend to explore the unit baryon number solutions 
the problem of regularizing the imaginary part plays a central role. Actually, a regularized 
imaginary part will not a priori yield an integer baryon number. 

The aim of this letter is to investigate the soliton solutions of the NJL model with a 
regularized imaginary part, i.e. a regularized baryon number. Unit baryon number will be 
enforced by adding an appropriate constraint to the energy functional. We will demonstrate 
that this constraint cannot be satisfied by considering only the chiral angle field. However, 
a unit regularized baryon number is attainable when additionally scalar degrees of freedom 
are allowed to vary in space. This feature is also related to the recently observed collapse 



of the NJL soliton |TI| . It is unstable against building up a narrow and infinitely high peak 
of the scalar field at the origin. The valence quark is thereby joining the Dirac sea such 
that infinitely many avoided crossings occur. This effectively corresponds to a situation 
where a level stemming from the positive part of the spectrum acquires an infinitely large 
negative energy eigenvalue. Without regularizing the baryon number such a configuration 
would correspond to B — 1. Obviously it should be clear that cutoff models like the NJL 
model cannot be trusted when infinitely large energies play a significant role. On the other 
hand, the regularized baryon number tends to zero during the collapse, i.e. baryon number 
is leaking out of the soliton. In this letter we will demonstrate that this leaking is prevented 
and the collapse is avoided by fixing the regularized baryon number to unity. Then stable 
solutions exist with scalar and pseudoscalar mesons includedP Furthermore we will discuss 
the influence of fixing the regularized baryon number on several baryon properties. 



The model 

The starting point of our calculations is the two-flavor NJL action Anjl which, after 
bosonization |fTH| , may be expressed as the sum Anjl = Af + A m of a fermion determinant 
and a purely mesonic part 



*The soliton with scalar and pseudoscalar mesons is known to be also stabilized by adding a four meson 
interaction jlj]. 
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A F = Tr log(ip) = Tr log (i$ - (P R M + P L M^j), 

An = J d A x ^~tr(M^M -m (M + M^) +m 2 Q ) S j . (1) 

Here, Pr^l are the projection operators on the right- and left-handed quark fields, respec- 
tively, and m = diag(m u ,m d ) denotes the current quark mass matrix. We will restrict 
ourselves to the two flavor case and assume isospin symmetry: m" = m d = m. The coupling 
constant g will be determined from meson properties. The complex field M. = S + iP de- 
scribes the scalar and pseudoscalar meson fields which can be parametrized by M. = $[/, $ 
and U being hermitian and unitary matrices, respectively. We will refer to $ as the chiral 
radius while the chiral angle G is introduced via U = exp(zG). 

The quark determinant Af diverges and must therefore be regularized. For the regu- 
larization procedure it is necessary to continue to Euclidean space. This yields a complex 
Euclidean action and we consider its real (Ar) and imaginary (Ai) parts separately: 

A R = ^Trlog(4|fe) 

M = ^Trlog^lr 1 ^)- (2) 
To keep the 0(4)-invariance, we use Schwinger's proper time regularization For the 



real part this prescription consists of replacing the logarithm by a parameter integral 

A R - -~ r -Trexp (-sp\p E ) , (3) 

2 Jl/A 2 s \ J 

which for A — ► oo reproduces the logarithm up to an irrelevant constant. 

Corresponding to the action (|l|) the energy of the static soliton E so1 = E s + E m splits 
into a fermionic and a mesonic part, where E f = E val + E vac contains valence and vacuum 
parts. To calculate W we express the Euclidean Dirac operator in its Hamiltonian form: 

ip E = P(-d T -h), (4) 

wherein r denotes the Euclidean time coordinate. Substituting the hedgehog ansatz U (r) = 
exp {it ■ rB(r)} for the chiral field and a radial function $ = $(r) for the chiral radius the 
static Hamiltonian reads: 

h = a ■ p + /3$ (cos 6 + ir ■ tj 5 sin 0) . (5) 

Denoting the associated eigenvalues by e M the fermionic contribution to the static energy is 
given by [110]: 



EmC t > $ ] = J™ ds s " 3/2 E H {-<) ~ ex P h« 



) 2 



N r 1 

/47T Jl/A* 

E™ l [Q,$] = N c J2v,K\ (6) 

where ?7 M denote the (anti-) quark occupation numbers and e° are the eigenvalues of h for 
the mesonic vacuum A4 = (A4) = Ml with M being the constituent quark mass. Note that 
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the vacuum energy has been subtracted in . The only free parameter is M since A is fixed 
by fitting the pion decay constant /„. = 93MeV. The gap equation relates the current quark 
mass m to the constituent quark mass M. The coupling constant g is eliminated by fitting 
the pion mass m n = 135MeV: g = mM/m 2 f 2 [jnj. Thus the mesonic part of the soliton 
energy may be expressed as 



E m [9, $] 



7T« 7T 

mM 



dr r z 



{$ 2 



M 



2m($(r) cos9(r) — M)| 



(7) 



Again the contribution of the vacuum configuration has been subtracted. 

Until now, we did not take into account regularization of Aj. Since this part is finite, the 
question arises whether it has to be regularized or not. Here we argue that regularization 
is necessary in order to have a consistent appearance of the one particle eigenenergies 
16[ O, IB, [18|. Ai does not contribute to the soliton mass, however, its regularization has 



drastic consequences when we regard the sea contribution to baryon number 



B vac = lim ^Tr \h (-d 2 + h 2 )- 1 ) 



where T is the Euclidean time interval under consideration. We introduce regularization of 
the baryon number by replacing 



'-dl + h 



2\-l 



1/A 2 



(9) 



which again is an identity for A — > oo. The regularized baryon number B\ then reads: 



wherein = — |erfc ( ^ ) denote the vacuum occupation numbers of the single quark 
orbits fi in the proper-time regularization 0. In the limit A — > oo, Ba is obviously 
integer, however, this is no longer the case for finite A. The main purpose of this paper is to 
investigate self-consistent solitons constrained to the baryon number B A = 1. Accordingly, 
the energy functional is replaced by 



E = E val + E vac + E m + M 



l) 2 - -a\ 2 
' 2 



'11' 



The constituent quark mass M is included such that the Lagrange multiplier A and the 
auxiliary parameter a are dimensionless. Varying the total energy functional E with respect 
to A yields 



A 



{Bp 



1) 



(12) 



The limit a — > has to be assumed in order to achieve B\ 
to the meson fields provides the equations of motion 



1. Variation of E with respect 



P(r) cos 9(r) 
$(r) 



sin 9(r) 



mN c M 

m cos 9(r) — -^-p i^v) cos 9(r) + P(r) sin 9(r)] (13) 
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where we have denoted 



S(r) = £ [(A/; + VlM ) sign(e M ) + C„] I dQ 

P( r ) = [{K + %) sign( eAt ) + C„] f dQ % [ir ■ r6(r) 75 ] (14) 

The derivative of the constraint with respect to the one particle eigenenergy is proportional 
to 
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_ 2MX(B A -1) f 6 
C " ~ iV c 0FA eXp \ A* J" (i5j 

Numerical Results 

The self-consistent soliton is constructed by iterating the constraint ([12]) and the equa- 
tions of motion ( |T5D together with the corresponding eigenvalue problem associated with the 
Hamiltonian (^). The latter is solved by discretizing the momentum eigenvalues of the free 
Hamiltonian. This is achieved by putting the system into a large spherical box of radius D. 
Our calculations are performed with a constituent quark mass M = 400 MeV. We consider 
this to be sufficient in order to demonstrate that the soliton is stabilized by constraining the 
regularized baryon number to unity. 

Restricting the meson profiles to the chiral circle and omitting the constraint ( pi]) we 
find B\ = 0.966 for the regularized baryon number. This number cannot be augmented 
considerably by including flllD and staying on the chiral circle. I.e. we do not find stable 
solutions in the limit a — > when the chiral angle is taken to be the only space dependent 
meson field. Allowing, in addition, the chiral radius to be space dependent as well, we 
face the well-known problem that the soliton of scalar and pseudoscalar fields collapses. As 
discussed in the introduction this collapse is associated with the leakage of baryon number 



B\ — > and as we will see, including the constraint flTT| ) stabilizes the soliton. Indeed we find 
stable solutions with scalar and pseudoscalar fields when the regularized baryon number is 
constrained to unity. As can be observed from table 1, B\ = 1 cannot exactly be fulfilled for 
finite D and we conjecture that B\ = 1 is only attainable in the continuum limit, D — > oo. 

In table 2 we display the energies of the self-consistent soliton solutions for various 
strengths of the constraint. These strengths are given in terms of the auxiliary parame- 
ter a and can be transferred to different regularized baryon numbers B\. One sees that the 
valence quark eigenenergy increases as B^ tends to unity. The resulting meson profiles are 
still localized as is obvious from figs. 1 and 2 where we display the meson profiles O(r) and 
$(r) corresponding to the regularized baryon numbers B^ of table 2. The chiral angle O(r) 
exhibits a strong squeezing at the origin with growing baryon number, while for r > 1 fm all 
profile functions show the same spatial dependence. Thus, the constraint effects the meson 
profiles only at small r. This can also be observed from the profile function $(r) which 
almost vanishes at the origin, reflecting a local partial restoration of chiral symmetry. 

In order to further investigate whether the solutions are localized objects we have con- 
sidered the baryon densities (cf. fig. 3). Here, the sea quark part of the baryon density is 
of particular interest because it almost vanishes, whereas the energy of the Dirac sea turns 
out to be sizable. In order to understand this result, we split in fig. 4 the sea quark density 
into two parts stemming from the intrinsic positive and negative parity eigenf unctions of (^), 
respectively. It can clearly be seen that the contributions of these parts to the baryon density 
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Fig. 1. The chiral angle @(r) for the non- linear NJL soliton (short dashed line) and different strengths of 
the constraint. 
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Fig. 2. Same as fig. 1 for the chiral radius $(r). 
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Table 1. The maximal baryon number obtained for different box radii D. 



D/fm 


B A 


4 


0.9970 


5 


0.9976 


6 


0.9977 


7 


0.9979 


8 


0.9980 



Table 2. The soliton energy E for constituent mass M =400 MeV with its valence quark, sea and mesonic 
contributions E val , E vac and E m for different baryon numbers Ba- The first line refers to the non-linear 
case. The small deviation of E from E val + E vac + E m is due to the constraint in (11). All energies are 
given in MeV. 



Ba 








E 


0.9662 


632 


572 


34 


1238 


0.9932 


914 


753 


-387 


1285 


0.9939 


925 


755 


-398 


1286 


0.9946 


941 


758 


-426 


1287 


0.9965 


996 


770 


-475 


1293 


0.9973 


1030 


776 


-510 


1294 


0.9978 


1043 


780 


-529 


1296 



cancel. For the static energy, however, they sum up coherently resulting in a relatively large 
value. 

Let us next turn to the discussion of a few static properties resulting from our soliton 
solution. In table 3 we display the isoscalar radius (rf^) 1 / 2 which is associated to the 
baryon charge distribution, the axial charge of the nucleon, g^, as well as the moment 
of inertia a 2 for collective rotations in isospace. The latter measures the A-nucleon mass 
difference M& — M n = 3/2a 2 . We should mention that the analytical form of these quantities 
in terms of sums involving matrix elements of the eigenstates of (|]) do not differ from the 
original treatment of the pure pseudoscalar case. We thus may employ the relevant formulas 



from refs. |T(], [19j . In agreement with the results for the static energy we observe for these 
quantities a dominating valence quark contribution as Ba tends towards unity. Considering 
e.g. the moment of inertia, a 2 , the Dirac sea contributions reduces from 23% to 4% when 



/ 2 \ 1 /2 2 
Table 3. Baryon radius (?"j_q) ' , axial charge gA and the moment of inertia Ct for different baryon num- 
bers. The first line denotes the data from the unconstrained soliton on the chiral circle ($ = M). 



Ba 


(rLo) 1/2 /f- 


aV(GeV)- 1 


9A 


0.9662 


0.77 


5.52 


0.72 


0.9933 


0.93 


9.11 


0.83 


0.9946 


0.96 


9.87 


0.84 


0.9965 


1.03 


11.94 


0.86 


0.9973 


1.08 


13.54 


0.87 


0.9978 


1.11 


14.62 


0.89 
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Fig. 3. The baryon density p[r) (full line) and its valence (short dashed line) and sea quark contribution 
(long dashed line) for the soliton with B\ = 0.9978. 
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Fig. 4. The sea quark contribution to the baryon density (full line) and its contributions from states with 
intrinsic positive parity (short dashed line) and negative intrinsic parity (long dashed line) for the soliton 
with B A = 0.9978. 
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we go from the pure pseudoscalar to B\ = 0.9978. This scheme may easily be understood 
by remarking that the quantities under consideration are rather sensitive to the large r- 
behavior of the fields: The valence quarks get enhanced of large r while the meson fields 
are squeezed at the origin. It is remarkable that the inclusions of the constraint improves 
our predictions for (r| =0 ) 1//2 and g& while a 2 ~ 15 GeV -1 is far off the empirical value 
a 2 = 3/2(M A - M n ) = 5 GeV" 1 . 

Conclusions 

The main result of this letter is that the collapse of the NJL soliton is prevented when 
constraining the regularized baryon number to unity. This constraint can only be fulfilled 
when the scalar meson field is allowed to be space dependent. Then it displays a partial 
restoration of chiral symmetry close to the origin. Despite of the fact that the valence quark 
energy eigenvalue approaches the one without soliton the valence quark wave function is still 
quite well localized, and the Dirac sea is still polarized. The vacuum contribution to the 
baryon number vanishes, however, the vacuum contribution to the soliton energy is sizeable. 
Whereas some static properties of the baryons, the isoscalar mean square radius and the 
axial coupling g^, are improved when considering the B\ = 1 soliton, the moment of inertia 
is too large and therefore the predicted nucleon-A splitting becomes too small. 

We have seen that constraining the regularized baryon number leads to a valence quark 
dominated picture of the soliton. It would therefore be interesting to include the constraint 
also in the case when vector and axial vector fields are present. As is well known, these 
models support the Skyrmion picture of the baryon, i.e. the valence quark joins the Dirac 
sea ||. First results indeed show that the inclusion of the constraint leads again to an 
increased valence quark eigenenergy which may even be positive. 

In this letter we used the proper time regularization which already "damps" the contri- 
bution of low lying levels to physical quantities, especially to the baryon number. From this 
point of view it is easily understandable that levels are pushed "out of the gap" . The question 
remains whether a regularization procedure which does only influence high-lying levels will 
lead to a different qualitative behavior. This idea can probably be tested using dimensional 
regularization which is known to even enhance contributions from low-lying levels EH . 
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